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A search for single top quark production is performed in the full e^p data sample 
t^J- ■ collected by the HI experiment at HERA, corresponding to an integrated luminosity of 

474 pb _1 . Decays of top quarks into a b quark and a W boson with subsequent leptonic or 
hadronic decay of the W are investigated. A multivariate analysis is performed to discrim- 
inate top quark production from Standard Model background processes. An upper limit 
on the top quark production cross section via flavour changing neutral current processes 
a(ep — > etX) < 0.25 pb is established at 95% CL. Limits on the anomalous coupling n tui 
are derived. 
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1 Introduction 



Top quarks are of particular interest in searches for new physics because their mass is close 
to the electroweak scale. In e ± p collisions at HERA the production of single top quarks is 
kinematically possible due to the large centre of mass energy of up to y/s = 319 GeV. Within the 
Standard Model (SM) the production of top quarks at HERA is however strongly suppressed. 
Therefore the observation of single top quark production would be a clear indication of new 
physics. In several extensions of the SM the top quark is predicted to undergo flavour changing 
neutral current (FCNC) interactions, which may lead to a sizeable top quark production cross 
section at HERA [1,2]. 

A search for single top quark production is performed using the full e ± p data sample col- 
lected by the HI experiment at HERA. The data correspond to an integrated luminosity of 
474 pb -1 , including 36 pb 1 of data taken at y^s = 301 GeV. The analysis is inspired by and 
supersedes the previous HI search for single top quark production based on an integrated lumi- 
nosity of 118 pb -1 [3]. Anomalous single top quark production also has been investigated by 
the ZEUS collaboration [4] and by the LEP experiments [5]. 

In the present analysis, single top quark production is detected via the decay of the top quark 
t — > bW. In the case of leptonic decays of the W boson, W — > Iv, the signature is a charged 
lepton (electron or muon) and missing transverse momentum, accompanied by a hadronic final 
state X with a high transverse momentum (P T ) due to the b quark. Events of this topology 
have been studied in a recently published analysis of events with isolated leptons and missing 
transverse momentum by HI [6]. In the case of the hadronic W boson decay channel, W — > qq, 
the signature of single top quark production consists of three high P T jets with an invariant mass 
compatible with the top quark mass. A multivariate discriminant based on a neural network is 
used to differentiate top quark production from SM background. A possible contribution from 
a top quark signal is extracted using the method of fractional event counting [7]. 

2 Single Top Quark Production 

The dominant process of SM single top quark production at HERA is via charged current (CC) 
production ep — > vtbX. The cross section for this process has been estimated [8] to be 0(1) fb 
and is thus not observable with the available integrated luminosity. Anomalous single top quark 
production within FCNC models, where the coupling of a top quark with an up-type quark U 
via a photon is described by a coupling K t u-y, is illustrated in figured] This process, as well as 
FCNC top quark interactions involving vector couplings to a Z boson V t uz, can be described 
by an effective Lagrangian [9]: 

c fcnc = J2 e ^n tUy ia, u A^U + —?—V 1U z*y l >UZ'> + h.c. , 
11 2A 2cosflw 

U=u,c 

where = (if 2) [7^,7"], 9 W is the Weinberg angle, e and g are the couplings to the elec- 
tromagnetic and weak gauge groups with £7(1) and 577(2) symmetries, respectively, ejj is the 
electric charge of up-type quarks, A^ u is the field strength tensor of the photon, is the vector 
field of the Z boson and A is a scale parameter. 
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Figure 1: Anomalous single top quark production in ep collisions via a FCNC coupling K tul . 

As the top quark mass is comparable to the ep centre of mass energy at HERA, the inter- 
acting parton in the proton must be at high Bj0rken-x. Contributions from the charm quark are 
therefore neglected (n tn = V tcZ = 0) since the c quark density in the proton is low at high 
Bj0rken-x. Similarly, the production of anti-top quarks is neglected, as this process involves 
sea anti-quarks in the initial state. 

The simulation of an anomalous single top quark signal is done using the event generator 
ANOTOP [3], which uses the leading order (LO) matrix elements of the complete e + q — > 
e + t—^e + b + W^e + b + f + f' process as obtained from the CompHEP program [10]. 
ANOTOP is used to calculate the production cross section and to study decays of top quarks in 
the HI detector. Only top quark decays t — > bW are considered as suggested by the strict limits 
on other possible top quark decays [11, 12]. 

The anomalous single top quark production cross section can then be parametrised as: 

a(ep -> etX, y/s) = c 7 • «4 7 + c z ■ V t 2 uZ + c jZ • «Wy • V tuZ . (1) 

The coefficients c 7 and c z are determined with ANOTOP at ^Js = 319 GeV. In this analysis 
m top is set to 175 GeV and by convention the scale parameter A is fixed to m top . If A = m top 
is reduced to 170 GeV, the cross section increases by 25%, mainly due to an increased phase 
space for the production of top quarks. This mass range encompasses the current top quark mass 
determination of CDF and D0 m top = 172.4 ± 1.2 GeV [13] and corresponds to a coefficient 
range from c 7 = 7.53 pb to 9.41 pb and c z = 0.26 pb to 0.32 pb. The values for c 7 include next 
to leading order (NLO) corrections [14], which increase the LO cross section by 17%. Including 
these NLO corrections results in an uncertainty related to the choice of the renormalisation 
and factorisation scales of about 5%. Since the contribution of the Z boson is small and no 
competitive sensitivity is expected, Z-exchange is neglected in this analysis and only the K tul 
coupling is considered. The interference term with coefficient c lZ contributes less than 1% of 
the total cross section and is therefore also neglected. 
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3 Standard Model Background Processes 



The contributions from SM processes to the background in the leptonic decay channels were 
studied extensively in the analysis of events with isolated leptons and missing transverse mo- 
mentum [6] and are only briefly described here. The main SM background process is the pro- 
duction of single W bosons. Further SM background processes may mimic the top quark sig- 
nature through misidentification or mismeasurement. In the hadronic decay channel the SM 
background prediction is dominated by multi-jet events produced in photoproduction or neutral 
current (NC) deep-inelastic scattering (DIS). 

Production of single W bosons is modelled by the EPVEC generator [15]. NLO QCD cor- 
rections are taken into account by reweighting the generated events [16, 17]. CC DIS events are 
simulated using the DJANGO [18] generator, which includes first order leptonic QED radiative 
corrections based on HERACLES [19]. The production of two or more jets in DJANGO is ac- 
counted for using the colour-dipole model [20]. The RAPGAP [21] event generator, which im- 
plements the Born, QCD Compton and Boson Gluon Fusion matrix elements, is used to model 
NC DIS events. Direct and resolved photoproduction of jets and prompt photon production 
are simulated using the PYTHIA [22] event generator. The simulation is based on Born level 
hard scattering matrix elements with radiative QED corrections. In RAPGAP and PYTHIA, 
jet production from higher order QCD radiation is simulated using leading logarithmic parton 
showers while hadronisation is modelled with Lund string fragmentation. The leading order 
MC prediction from NC DIS and photoproduction processes with two or more high transverse 
momentum jets is scaled by a factor of 1.2, to account for missing higher order QCD contribu- 
tions in the MC generators [23,24]. Further small contributions to the SM background originate 
from lepton pair production, simulated using the GRAPE generator [25], and from photon wide 
angle bremsstrahlung, simulated in the WABGEN framework [26]. 

Generated signal and background events are passed through the GEANT [27] based simu- 
lation of the HI detector, which takes into account the running conditions of the different data 
taking periods, and are reconstructed and analysed using the same program chain as for the data. 



4 Experimental Conditions 

A detailed description of the HI experiment can be found in [28]. Only the detector compo- 
nents relevant to the present analysis are briefly described here. The origin of the HI coordi- 
nate system is the nominal ep interaction point, with the direction of the proton beam defin- 
ing the positive z-axis (forward region). Transverse momentum is measured in the rry-plane. 
The pseudorapidity 77 is related to the polar angle 9 by 77 = — In tan(0/2). Tracking is pro- 
vided by the central (20° < 6 < 160°) and forward (7° < 6 < 25°) tracking detectors. They 
are used to measure charged particle trajectories, to reconstruct the interaction vertex and to 
complement the measurement of hadronic energies. The tracking detectors are surrounded 
by a finely segmented Liquid Argon (LAr) calorimeter [29] that covers the polar angle range 
4° < 6 < 154°. Electromagnetic shower energies are measured with a precision of a(E)/E = 
12%/y/E/GeV © 1% and hadronic energies with a(E) /E = 50%/^/E/GeV © 2%, as mea- 
sured in test beams [30, 31]. In the backward region, energy measurements are provided by a 
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lead/scintillating-fibre (SpaCal) calorimeteiQ [32] covering the range 155° < 9 < 178°. The 
LAr calorimeter and inner tracking detectors are enclosed in a super-conducting magnetic coil 
with a field strength of 1.16 T. From the curvature of charged particle trajectories in the mag- 
netic field, the central tracking system provides transverse momentum measurements with a 
resolution of op T / Pt = 0.005Pr/GeV © 0.015 [33]. The return yoke of the magnetic coil 
is the outermost part of the detector and is equipped with streamer tubes forming the central 
muon detector (4° < 9 < 171°). In the forward region of the detector (3° < 9 < 17°) a set of 
drift chambers detects muons and measures their momenta using an iron toroidal magnet. The 
luminosity is determined from the rate of the Bethe-Heitler process ep — * ep'j, measured in a 
photon detector located close to the beam pipe at z = —103 m, in the backward direction. 

The LAr calorimeter provides the main trigger for events in this analysis. The trigger effi- 
ciency is almost 100% for events with an electron with an energy above 10 GeV. Events with 
muons are triggered by an imbalance in transverse momentum measured in the calorimeter 
P£ al °. The trigger efficiency is ~ 60% for P£ al ° = 12 GeV, rising to ~ 98% for P£ al ° > 25 GeV. 
Events with jets have a trigger efficiency of nearly 100% for events with Pf ct > 25 GeV [34]. 

In order to remove background events induced by cosmic rays and other non-ep sources, the 
event vertex is required to be reconstructed within ±35 cm in z of the average interaction point. 
In addition, topological filters and timing vetoes are applied. 

5 Data Analysis 

5.1 Particle Identification and Event Reconstruction 

Particle identification and hadronic final state reconstruction are described in detail in [6] and 
are only briefly summarised here. Electrons and photons are characterised [35] by compact 
and isolated electromagnetic clusters in the LAr calorimeter or SpaCal. Muon identification is 
based on a track in the inner tracking systems associated with signals in the muon detectors [36, 
37]. Calorimeter energy deposits and tracks not previously identified as electron, photon or 
muon candidates are used to form combined cluster-track objects, from which the hadronic final 
state is reconstructed [34,38]. Jets are reconstructed from these combined cluster-track objects 
using an inclusive kx algorithm [39, 40] with a minimum transverse momentum of 4 GeV. 
The missing transverse momentum P™ 1SS in the event is derived from all detected particles and 
energy deposits in the event. 

Strict isolation criteria are applied to electron and muon candidates. For electrons a track 
with a distance of closest approach (DCA) to the cluster of less than 12 cm is required to reject 
photons. The calorimetric energy measured within a distance in the pseudorapidity-azimuth 
(rj — (f)) plane D = ^/ Arj 2 + A0 2 < 0.5 around the electron cluster is required to be below 
3% of its energy. Furthermore, electrons are required to be isolated from jets by a distance 
D(e; jet) > 1.0 to any jet axis, and in the region 9 > 45° by a distance D(e; track) > 0.5 from 
any track. A muon candidate may have no more than 5 GeV energy deposited in a cylinder, 
centred on the muon track direction, of radius 25 cm and 50 cm in the electromagnetic and 
hadronic sections of the LAr calorimeter, respectively. Muon candidates are required to be 
separated from any jet by a distance D(fi; jet) > 1.0 and from any track by D(n; track) > 0.5. 

'This device was installed in 1995, replacing a lead-scintillator "sandwich" calorimeter [28]. 
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5.2 Systematic Uncertainties 



The following experimental systematic uncertainties are considered: 

• The uncertainty on the electromagnetic energy scale varies depending on the polar angle 
from 0.7% to 2%. The polar angle measurement uncertainty is 3 mrad for electromagnetic 
clusters. 

• The scale uncertainty on the transverse momentum of high P T muons is 2.5%. The un- 
certainty on the muon polar angle measurement is 3 mrad. 

• The hadronic energy scale is known within 1.5% [41]. The uncertainty on the jet polar 
angle measurement is 10 mrad [34]. 

• The uncertainty on the trigger efficiency is 3%. 

• The luminosity measurement has an uncertainty of 3%. 

The effect of the above systematic uncertainties on the MC expectation is determined by varying 
the experimental quantities by ±1 standard deviation in the MC samples and propagating these 
variations through the whole analysis. 

Additional model uncertainties are attributed to the SM MC generators described in sec- 
tion[3] A theoretical uncertainty of 15% is used for the predicted contributions from EPVEC [16]. 
In the electron and muon channels the CC background contribution, which is modelled using 
DJANGO, is attributed a systematic error of 50% [6]. The contributions from background pro- 
cesses modelled using RAPGAP, PYTHIA, GRAPE and WABGEN are attributed 30% model 
uncertainties [6]. In the hadronic channel, the background processes with multi-jet final states 
modelled using RAPGAP and PYTHIA are attributed a 10% model uncertainty determined 
from the comparison to data in an extended phase space. 

5.3 Electron and Muon Channels 

The search for single top quark production in the electron and muon channels is based on the 
selection described in [6]. Isolated electrons and muons with a transverse momentum P^ > 
10 GeV in the polar angle range 5° < 0£ < 140° are selected in events with a missing transverse 
momentum P™ 1SS > 12 GeV. In the electron channel, 39 events are found, compared to a 
SM prediction of 43.1 ± 6.0. In the muon channel, 14 events are observed, compared to a 
SM prediction of 11.0 ± 1.8. To estimate a potential top contribution to this sample, a top 
quark candidate is reconstructed from its decay products (lepton £, neutrino v and b quark), 
and the compatibility with single top quark production via FCNC is tested using a multivariate 
discriminant method. 

The neutrino four-vector P v is reconstructed using the transverse and longitudinal momen- 
tum balance of the event. The transverse momentum of the neutrino is reconstructed by assum- 
ing: 

p/y prniss 

1 rp — T 
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If the scattered electron is detected in addition to the isolated lepton, the four-vector of the neu- 
trino can be fully reconstructed by exploiting the energy and longitudinal momentum balance: 

J2 - p l) + ( EU - p z) = 2E e = 55.2 GeV, (2) 

i 

where the sum runs over all detected particles, P z is the momentum along the proton beam 
axis and E® is the electron beam energy. For events with more than one isolated electron, the 
electron with the lower transverse momentum is assumed to be the scattered electron [42]. If 
the scattered electron is not detected, the constraint M lv ~ Mw = 80.42 GeV is applied. In 
the case that two physical solutions are obtained for [E v — P"), the more probable solution 
according to the ANOTOP simulation is chosen [42]. A small fraction of events is removed in 
the case when the reconstruction algorithm finds no physical solution. In the electron channel, 
38 events are selected after this neutrino reconstruction, while 39.7 ± 5.6 are expected from 
the SM. In the muon channel, 13 events are selected, while 10.7 ± 1.6 are expected from the 
SM. The selection efficiency for generated ANOTOP events at this stage is 49% (44%) for the 
electron (muon) channel. The four- vector of the 6-jet candidate is defined as the four-vector of 
the hadronic final state. The four- vector of the top quark candidate is reconstructed as the sum 
of the four- vectors of the isolated lepton, the neutrino and the 6-jet candidate. This includes top 
quark decays where the W boson decays via W — > r — > e(/i). Figure [2] shows the reconstructed 
top mass in the combined electron and muon channels. The data are in overall agreement with 
the SM prediction, while a slight excess of data events is observed in the top quark mass range. 

In addition to the kinematic reconstruction of the top quark decay the lepton charge is also 
exploited, as the decay chain t — > bW + — > iub produces only positively charged leptons. Well 
measured negatively charged leptons are rejected by requiring qe ■ > —1.0 where qi is the 
charge of the lepton, \k\ is the curvature of the track associated to the lepton and 5k is the 
error on the curvature. This requirement is only applied in the central region where the charge 
determination by the tracking detectors is reliable. Table \T\ summarises the event yield in the 
electron and muon channels for the resulting top preselection. In the electron channel 30 events 
are selected, while 31.5 ± 4.0 are expected from the SM. In the muon channel, eight events are 
selected, while 8.0 ± 0.9 are expected from the SM. 

The following observables are then used to further discriminate single top quark production 
via FCNC against SM background: 

• PJp, the transverse momentum of the 6-jet candidate; 

• M Eub , the invariant mass of the reconstructed top quark; 

• 9^, , the W decay angle calculated as the angle between the lepton momentum in the W 
rest frame and the W direction in the top quark rest frame. 

Distributions of these observables are shown in figures [3] (a-c) and|4](a-c) for the electron and 
muon channels, respectively. Good overall agreement of the data with the SM expectation is 
seen in all distributions. For the discrimination of top and background events, a multivariate 
discriminator based on a multilayer perceptron (MLP) neural network was trained on the sig- 
nal and background simulations. Figures [3] (d) and@](d) show the MLP discriminator output 
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distributions for the electron and muon channels, respectively. A good agreement is observed 
between data and SM expectation, with some events visible in the signal region. 

Results from the multivariate analysis are cross checked with a cut-based top selection re- 
quiring for the top preselected events in addition Pj, > 30 GeV and M iub > 140 GeV. The 
resulting event yields are also shown in tabled] In this selection five events are selected in the 
electron channel, while 3.2 ± 0.4 are expected from the SM. In the muon channel, four events 
are selected, while 2.1 ± 0.3 are expected from the SM. This sample includes all five top-like 
events found in the previous analysis [3]. If the charge requirement is not applied, the corre- 
sponding event yields are seven selected events for an expectation of 4.1 ± 0.7 in the electron 
channel, and six events for 2.8 ± 0.4 expected in the muon channel. 

5.4 Hadronic Channel 

The hadronic decay of the W boson from the top decay t — > bW — > bqq' leads to events with 
at least three jets with high transverse momenta. SM background arises mainly from multi-jet 
events in photoproduction or NC DIS. A top preselection is defined by selecting events with 
at least three jets in the pseudorapidity range —0.5 < rf et < 2.5. The jets are ordered by the 
magnitude of their transverse momenta and only events with P^ ctl > 40 GeV, P^ t2 > 30 GeV 
and Pf t3 > 15 GeV are selected. A cut on the scalar sum of the jet transverse momenta 
> HO GeV is also applied. In addition, one of the jet pairings must yield an invariant 
mass between 65 GeV and 95 GeV, corresponding to a window around the nominal W mass 
with a width of twice the mass resolution obtained for hadronic W decays [42]. The remaining 
jet that is not used to form the W candidate is considered to originate from the b quark and 
is required to have a minimum P T > 25 GeV. The yield for this top preselection is given in 
table CD A total of 404 events are selected, compared to a SM prediction of 388 ± 32 events. 
The selection efficiency for generated ANOTOP events is 42% at this stage. Analogously to 
the leptonic channel observables for the discrimination of the top quark signal from the QCD 
background are chosen: 

• P|,, the transverse momentum of the 6-jet candidate; 

• M jets , the invariant mass of all jets in the event corresponding to the mass of the top quark 
for signal events; 

• Oyy* , the W decay angle defined as the angle in the W rest frame between the lowest Pt 
jet of the two jets associated to the W decay and the W direction in the top quark rest 
frame. 

Distributions of these three observables are shown in figure \5\ (a-c), compared to the SM expec- 
tation and the simulated top signal. Good agreement between the data and the SM simulation 
is seen for all three distributions. For the discrimination of top and background events, an MLP 
discriminator was trained on the signal and background simulations. Figure |5] (d) shows the 
discriminator output distribution for the signal and background simulations. Also shown are the 
data events as classified by the discriminator and a good agreement is observed between data 
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and SM expectation. This distribution is used to extract a possible signal, as described in the 
next section. 

A cut-based top selection is also performed in the hadronic channel. The transverse mo- 
mentum of the 6-jet candidate has to fulfill P% > 40 GeV and the reconstructed top quark mass 
150 < Af jcts < 210 GeV. The number of candidate events selected is 128, compared with 
123 ± 13 events expected from SM processes. 

6 Results 

The discriminant distributions observed in all three channels are in correspondence with the 
SM expectation as shown in figures [3] (d),@](d) and[5](d). Since no significant indication for 
the production of single top quarks is observed, upper limits on the cross section are derived 
using the method of fractional event counting which takes efficiencies, statistical and system- 
atic uncertainties into account [7]. This method is cross-checked using a modified frequentist 
approach [43] which is found to yield similar results. Also expected cross section limits are 
determined by assuming that a top quark signal is a fluctuation of Poisson distributed SM back- 
ground. An ensemble of toy experiments with this hypothesis is conducted using the back- 
ground prediction and calculating a cross section limit at 95% CL for each toy experiment. The 
mean of the resulting distribution is reported as the expected cross section limit at 95% CL. 

For all channels combined an upper bound on the cross section for single top quark produc- 
tion via FCNC is obtained at the 95% CL: 

a(ep -> etX, y/s = 319 GeV) < 0.25 pb. 

This cross section limit is reported for y/s = 319 GeV, taking into account the ratio 0.70 of the 
predicted signal cross sections at y/s = 301 GeV and 319 GeV [14]. Table[2]lists the observed 
and expected cross section limits for the individual channels, the combined leptonic channels, 
and all channels combined. 

The top-like events observed in the data at high values of the discriminant lead to observed 
limits that are weaker than the expected limits. This difference corresponds to the slightly higher 
observed event yields from the cut-based top selection compared to the predictions. Extracting 
a production cross section from the discriminant distributions gives for all channels combined 
a(ep — > etX, y/s = 319 GeV) = 0.11 ± 0.07 pb. This value is compatible with zero within 
two standard deviations. 

The limits on the cross section are converted to limits at 95% CL on the anomalous FCNC 
coupling K tuJ using equation [TJ This results in an upper limit on K tul < 0.16 — 0.18 for the 
parameter range A = m top = 170 — 175 GeV, respectively. Figure [6] shows current limits 
from experiments in the n tul -V tu z plane. The limits on branching fractions of anomalous 
decays of top quarks obtained by the CDF experiment at Tevatron [11, 12] are converted to 
limits on couplings using the conventions described in section 2. The limit on the branching 
fraction B(t — > qZ) < 3.7% at 95% CL reported by CDF [11] corresponds to the strictest 
limit on the V tu z coupling. For ease of comparison, the HI limit on the coupling K tui is also 
converted to a limit on the branching fraction B(t — > try) < 0.64%, also shown in tablefU The 
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LEP experiments have also searched for anomalous single top quark production [5]. The limit 
obtained by the ZEUS experiment [4] using 130 pb' 1 of data is similar to the HI bound based 
on the full HERA data, mainly reflecting the slight excess of top-like events observed by HI, 
but not by ZEUS. The present analysis explores a domain not covered by other colliders. 



7 Summary 

A search for single top quark production via FCNC at HERA is performed using the full e ± p 
data sample, corresponding to an integrated luminosity of 474 pb _1 . The search is based on a 
sample of events with isolated leptons and missing transverse momentum and a sample of multi- 
jet events. No evidence for single top quark production is observed. Within FCNC models the 
top quark production cross section is extracted from multivariate discriminator distributions 
using the method of fractional event counting. An upper limit at 95% CL on the production 
cross section: 

a(ep -> etX, y/s = 319 GeV) < 0.25 pb 
is established. This result is translated into a limit on the anomalous coupling K tuj at 95% CL: 

K tvrt < 0.18 

if the scale of the new physics A = m top = 175 GeV. This corresponds to a limit on the 
branching ratio B(t — > wy) < 0.64%. 
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HI Search for Single Top Production at HERA (e±p, 474 pb _1 ) 
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2.1 ±0.3 
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Hadron Channel 


Data 


Standard Model 


Top Efficiency 
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404 


388 ± 32 


42% 


Cut-based Top Selection 


128 


123 ± 13 


33% 



Table 1: Observed and predicted numbers of events in the electron, muon and hadronic channels 
for the steps of the top preselection and in the cut-based top selection. The top signal efficiency 
is estimated using ANOTOP. The quoted errors take into account all statistical and systematic 
uncertainties. 
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HI Search for Single Top Production (e ± p, 474 pb x ) 


Channel 


Upper Limit at 95% CL 


a(ep -^tX,y/s = 319 GeV ) 
Observed [pb] Expected [pb] 




B(t -> wy) 
[%] 


Electron 
Muon 
Electron+Muon 
Hadronic 
Combined 


0.40 0.24 
0.30 0.22 
0.27 0.15 
0.42 0.27 
0.25 0.12 


0.21 - 0.23 
0.18-0.20 
0.17-0.19 
0.21 - 0.24 
0.16-0.18 


0.82 - 1.02 
0.61-0.76 
0.55-0.69 
0.86-1.07 
0.51 - 0.64 



Table 2: Observed and expected upper limits at 95% CL on the single top production cross 
section at ^J~s = 319 GeV in the electron, muon, hadronic and combined channels. Also shown 
are upper limits on the anomalous coupling n tuJ and on the branching fraction B(t — > wy). 
These limits are shown for A = m top = 170 — 175 GeV. 
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Figure 2: Distribution of the reconstructed top mass M in the electron and muon channels 
after v reconstruction but before the cut on the lepton charge. The data are shown as points, the 
total SM expectation as the open histogram with systematic and statistical uncertainties added in 
quadrature (shaded band). Also shown is the ANOTOP prediction with arbitrary normalisation 
(dashed histogram). 
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Search for Single Top Quark Production at HERA (474 pb" 1 ) 




Figure 3: Distributions of observables in the electron channel used to differentiate single top 
production from SM background processes in the top preselection. Shown are the reconstructed 
top mass M evb (a), the transverse momentum of the reconstructed 6-jet candidate P% (b), the W 
decay angle cos 9%y* (c), and the combination of these observables into a MLP-based discrimi- 
nant D (d). The data are shown as points, the total SM expectation as the open histogram with 
systematic and statistical uncertainties added in quadrature (shaded band). Also shown is the 
signal prediction with arbitrary normalisation (dashed histogram). 
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Search for Single Top Quark Production at HERA (474 pb" 1 ) 
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Figure 4: Distributions of observables in the muon channel used to differentiate single top pro- 
duction from SM background processes in the top preselection. Shown are the reconstructed 
top mass M^ ub (a), the transverse momentum of the reconstructed 6-jet candidate (b), the 
W decay angle cos Qyy* (c), and the combination of these observables into a MLP-based dis- 
criminant D (d). The data are shown as points, the total SM expectation as the open histogram 
with systematic and statistical uncertainties added in quadrature (shaded band). Also shown is 
the signal prediction with arbitrary normalisation (dashed histogram). 



20 



Search for Single Top Quark Production at HERA (474 pb" 1 ) 
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Figure 5: Distributions of observables in the hadronic channel used to differentiate single top 
production from SM background processes in the top preselection. Shown are the invariant 
mass of all jets M jets (a), the transverse momentum of the 6-jet candidate P|> (b), the W decay 
angle cos 6^, (c), and the combination of these observables into a MLP-based discriminant 
D (d). The data are shown as points, the total SM expectation as the open histogram with 
systematic and statistical uncertainties added in quadrature (shaded band). Also shown is the 
signal prediction with arbitrary normalisation (dashed histogram). 
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Limits on Single Top Quark Production via FCNC 




Figure 6: Upper limits at 95% CL on the anomalous K tui and V tu z couplings obtained as- 
suming A = m top = 175 GeV. Anomalous couplings of the charm quark are neglected: 
= Vtcz = 0. The domain excluded by HI is represented by the light shaded area. 
The dark shaded band shows the region additionally excluded if A = m top is varied to 
170 GeV. The hatched area corresponds to the domain excluded at the Tevatron by the CDF 
experiment [11, 12]. Also shown are limits obtained at LEP by the L3 experiment [5] and at 
HERA by the ZEUS experiment [4]. 
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